Mean-field potential calculations of high-pressure equation of state for 

shock-compressed BeO 
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A systematic study of the Hugoniot equation of state, phase transition, and the other thermody- 
namic properties including the Hugoniot temperature, the electronic and ionic heat capacities, and 
the Griineisen parameter for shock-compressed BeO, is presented by calculating the total free en- 
ergy. The method of calculations combines first-principles treatment for 0-K and finite-T electronic 
contribution and the mean-field-potential approach for the vibrational contribution of the lattice ion 
to the total energy. Our calculated Hugoniot shows good agreement with the experimental data. 
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As one member of the series of alkaline-earth ox- 
ide, Beryllium oxide (BeO) has many unique properties. 
First, it is the only one alkaline-earth oxide to crystal- 
lize in the wurtzite (WZ) rather than the rock-salt (RS) 
structure. Further, not only is BeO harder than the other 
alkaline-earth oxides but also it is among the hardest ma- 
terials known. It is a good insulator like other alkaline- 
earth oxide, but its heat conductivity is an order of mag- 
nitude higher, which makes it a technically promising 
ceramic 1 . Also, BeO has a high melting point. These 
interesting physical properties are related to characteris- 
tic features of interatomic bonding in BeO: for example, 
Compton scattering measurements revealed a significant 
component of the primarily ionic bonding in BeO 2 . 

An important issue for BeO is its pressure induced 
phase transition sequence, which has been investigated 
by a few experimental^ and theoretical studies^i^. 
This is motivated by the well-known dielectric theory of 
Phillips and Van Vechten 1 ^ that the structural phase sta- 
bilities of the A N B S ~ N compounds can be related to their 
ambient spectroscopic ionicities (fi). According to this 
theory, binary compounds with /j > 0.785 crystallize in 
a sixfold coordinated structure such as RS, while those 
with fi < 0.785 crystallize in a fourfold structure such 
as zincblende (ZB) or wurtzite (WZ). BeO has a Phillips 
ionicityiS of 0.602 and crystallize in the WZ structure 
under ambient conditions, in accord with the dielectric 
theory. Under pressure, the dielectric theory^ predicts 
that the tetragonal compounds with fi > 0.35 first trans- 
form to an ionic sixfold structure and then to some metal- 
lic phase, while those with /, < 0.35 transforms to the 
metallic phase directly. Given this understanding, how- 
ever, the pressure induced phase sequence for BeO is 
not yet fully understood. The early two ab initio cal- 
culations support this general prediction that a WZ— >RS 
transition occurs at 22 GPa£ (40 GPa from Ref. 4 ). On 
the other side, the latter calculations predicted the WZ 
first transformed into zinc blende (ZB) and then into RS. 
Herein, the work of Camp and Doren using soft non- 
local pseudopotential predicted the WZ-ZB-RS transi- 
tions at 74 and 137GPaii, while Park et al^, who used 
a first-principles soft nonlocal pseudopotential method 
within the generalized-gradient approximation (GGA), 



predicted the transition pressure to be 91 (WZ— >ZB) and 
147 GPa (ZB— >RS) respectively. This unusual transition 
sequence is attributed to the large charge asymmetry and 
the small bond length of BeOi£. By use of full-potential 
linear muffin-tin orbital (FPLMTO) method, Boettger 
and Wills^ have also predicted a WZ^ZB^RS phase 
transition sequence with the final transition at about 95 
GPa. They also predicted a theoretical Hugoniot, consis- 
tent with shock wave data up to 100 GPa. More recently, 
Cai et al*£, by systematically calculating the enthalpy 
barrier of the phase transition, have shown that only the 
WZ— >RS transition will occur with increasing pressure. 

Experimentally, Raman spectra showed no evidence for 
a phase transition up to 55 GPa in BeO 4 and Hugoniot 
data did not reveal any volume change for stresses up to 
100 GPa 4 . The most recent static high pressure x-ray 
diffraction experiment^ also showed no phase transition 
up to 126 GPa, but a new phase was found at 137 GPa. 
With increasing pressure, the WZ phase disappeared at 
175 GPa.ii 

In this paper we theoretically study the phase transfor- 
mation and Hugonoit of BeO under shock-wave dynamic 
compression, which is one of the most efficient ways to 
explore the thermodynamic properties of a material at 
ultrahigh pressures and temperatures. From its derived 
Hugoniot state one can deduce much useful information 
such as equation of state (EOS) at feasible temperature 
range. The theoretical methodology we employed in this 
paper is a combined ab initio electronic structure cal- 
culation and classical mean-field potential (MFP) ther- 
modynamic treatment 1 ^. Recent efforts have applied the 
MFP approach to a number of prototype metals 1 ^, in- 
dicating that both the calculated Hugoniot equation of 
state (at pressures up to ^1000 GPa and temperatures 
up to ^70 000 K) and room temperature isotherms are 
well described within the experimental uncertainties. We 
calculate the OK enthalpy of WZ, ZB, and RS BeO as a 
function of pressure. The results show that the enthalpy 
of WZ BeO is equal to that of ZB BeO at 84 GPa. When 
pressure is increased to be 105 GPa, the enthalpy for 
WZ and RS BeO approach to cross. Based on our first- 
principles total energy data, we give a MFP analysis of 
the phase transition, Hugoniot EOS and the other ther- 
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modynamic properties for shock-compressed BeO. Our 
calculated Gibbs free energies show that the WZ— »RS 
trnasition pressure is 103.8 GPa. Furthermore, our cal- 
culated Hugoniot is in good agreement with the experi- 
mental data. 

For a system with a given averaged atomic volume V 
and temperature T, the Helmholtz free-energy F(V, T) 
per atom can be written as 

F(V, T) = E C (V) + F lon {V, T) + F el (V, T), (1) 

where E c represents the 0-K total energy which is 
obtained from ab initio electronic structure calcula- 
tions, F e i is the free energy due to the thermal excita- 
tion of electrons, and Fi on is the ionic vibrational free 
energy which is evaluated from the partition function 
Zion — exp(— NF ion /kBT). Here N is the total num- 
ber of lattice ions. In the mean-field approximation, the 
classical Z ion is given byi£ 
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The essential of the MFP approach is that the mean-field 
potential g(r, V) is simply constructed in terms of E c as 
follows^ 

g(r, V) = \ [Ec(R + r) + E C (R - r) - 2E C {R)] , (3) 

where r represents the distance that the lattice ion 
deviates from its equilibrium position R. It should 
be mentioned that the well-known Dugdale-MacDonald 
expression 17 for the Griineisen parameter can be derived 
by expanding g(r,V) to order r 2 . Then, F ion can be 
formulated as 



F on (V,T) 



k L J\\\n r ^-+\nv f {V,T)). (4) 



with 



v f (V,T) = 4ttJ exp 



9(r,V) 



r 2 dr. 



(5) 



When the electron-phonon interaction and the mag- 
netic contribution are neglected, the electronic contribu- 
tion to the free energy is F e i = E e i — TS e i, where the 
bare electronic entropy S e i takes the formic 

S el (V, T) = -k B J n(e, V) [/ In / + (1 - /) ln(l - /)] de, 

(6) 

where n(e, V) is the electronic density of states (DOS) 
and / is the Fermi distribution. With respect to Eq. 
(6), the energy E e i due to electron excitations can be 
expressed as 



E el (V,T)= n(e,V)fede- n(e,V)ede, (7) 
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FIG. 1: (color online). Calculated total energies for 
WZ,ZB,and RS phase of BeO are plotted as a function of 
volume. 



where ep is the Fermi energy. Given the Helmholtz free- 
energy F(V, T), the other thermodynamic functions such 
as the entropy S = —(8F/dT)v, the internal energy E — 
F + TS, the pressure P = -{dF/dV) T , and the Gibbs 
free energy G = F + PV, can be readily calculated. 

The first task is to calculate the 0-K total energy 
E C (V) in a wide range of the volume. For this we em- 
ploy the full-potential linearized-augmented-plane-wave 
(FPLAPW) method 20 within the generalized gradient 
approximation (GGA) 21 . The total energy of BeO with 
WZ structure was calculated at 96 lattice constants rang- 
ing from 4.14 a.u. to 6.04 a.u. in steps of 0.02 a.u., while 
for the ZB structure 91 lattice constants ranging from 
6.22 a.u. to 8.02 a.u. in steps of 0.02 a.u. were con- 
sidered. Finally for RS BeO we sample 81 lattice con- 
stants ranging from 5.46 a.u. to 8.66 a.u. in steps of 
0.04 a.u.. As a reference, we mention that the calculated 
zero-pressure lattice constants for the three BeO struc- 
tures are a = 5.135 a.u. (c/a = 1.623, WZ), 7.232 a.u. 
(ZB), and 6.894 a.u. (RS), respectively. In all of three 
structures of BeO and for all the atomic volumes con- 
sidered, we use the constant value of the muffin-tin radii 
R mt of 1.2 a.u. for Be atom and 1.0 a.u. for O atom 
respectively. The plane-wave cutoff K cut is determined 
by R m t x K cut — 10.0. 4000 fc-points in the full Brilloum 
zone are used for reciprocal-space integerations. 

Figure 1 shows the calculated 0-K total energies as a 
function of volume for WZ, ZB, and RS phases of BeO. 
The equilibrium volumes (per BeO molecule) correspond- 
ing to the minima in total energy are 95.53 a.u. 3 (WZ), 
95.02 a.u. 3 (ZB), and 82.18 a.u. 3 (RS), respectively. 
The equilibrium properties, including lattice parameters, 
bulk modulus from the fit with Murnaghan curved 2 ., and 
the equilibrium elastic constants for the three crystal 
structures of BeO, are listed in Table I. For comparison, 
the available experimental data are also given in Table I. 
The other available theoretical data have been collected 



3 



_0 

=3 
O 



X 




-0.0005 



20 40 60 80 100 120 140 160 180 200 

P (GPa) 



FIG. 2: Calculated enthalpy difference between WZ and ZB 
type BeO at K as a function of pressure. 
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FIG. 3: (color online). Calculated enthalpy of WZ and RS 
type BeO at K as a function of pressure. 



in Ref . [2] . For WZ phase, the calculated equilibrium lat- 
tice constant and bulk modulus are 2.712 Aand 203 GPa, 
in good agreement with the experimental values of 2.698 
Aand 212 GPa 3 . Also one can see that the calculated 
values of elastic constants c^- for WZ structure agree as 
well with the measured data. 



TABLE I: Calculated and experimental structural parame- 
ters for WZ, ZB, and RS phases of BeO: equilibrium lattice 
constants a (A), bulk modulus Bo (GPa), and the elastic con- 
stants Cij (GPa). Note that the c/a ratio for the WZ has been 
fully optimized. 



Phase 


a 


c/a 


Bo 


Cll 


C12 


Cl3 


C33 


C44 


C66 


WZ 


2.712 


1.623 


203 


432 


120 


87 


463 


142 


156 


WZ (exp. a ) 


2.698 


1.622 


212 


460 


125 


82 


490 


145 


167 


ZB 


3.825 




201 


342 


148 






208 




RS 


3.648 




231 


298 


214 






294 





1 Reference 23. 



To obtain the most stable structure at finite pressure 
and temperature, the Gibbs free energy G — E+PV—TS 
should be used. For the experimental data, the last 
term is assumed small and therefore usually neglected by 
the previous theoretical studies^&i^. For comparison 
with those calculations, we also work with the enthalpy 
H = E + PV for each phase as a first step. The enthalpy 
of BeO with WZ, ZB, and RS structures are calculated, 
and the phase transition pressure is obtained from the 
enthalpy curve crossings. The difference in enthalpy be- 
tween WZ and ZB type BeO as a function of pressure is 
plotted in Fig. 2, which shows that the two structures are 
energetically very close with the enthalpy difference typ- 
ically of ^5 meV in a wide range of pressure. This result 
is not surprising because of the high similarity between 
the two structures. The local environment of any atom in 



either ZB or WZ is exactly the same as far as the second 
neighbors. One can see from Fig. 2 that the enthalpy of 
WZ BeO is lower than that of ZB phase up to 84 GPa of 
pressure, while the enthalpy of ZB BeO begins to be lower 
for P > 84 GPa. It thus seems that the zero tempera- 
ture WZ^ZB transition pressure is 84 GPa. This result 
is close to that (87 GPa) obtained by Cai et alJ£ based 
on pseudopotential method with generalized-gradient ap- 
proximation (GGA). From the point of view of enthalpy 
barrier, however, Cai et al. noticed that the WZ— >ZB 
phase transition cannot happen up to 200 GPa. Here 
we adopt this point and will not consider the WZ— >ZB 
transition in the following Hugoniot calculations. 

The enthalpy of BeO with WZ and RS structures are 
shown in Fig. 3 as a function of pressure. One can see 
that the enthalpy of the two structures increase almost 
linearly with increasing pressure. When the pressure is 
up to 105 GPa they become equal. With further increas- 
ing pressure the RS structure becomes more stable. Thus 
it is expected to undergo a WZ^RS phase transition 
around P=105 GPa. This value is identical to that ob- 
tained in Refii 3 - but prominently differs from the calcu- 
lated result of 147 GPa in RefJ^. This difference seems 
to be mainly caused by the different methods for finding 
the transition pressure. 

Next we investigate the phase transition properties of 
BeO under shock-wave compress. In this case, we fully 
take into account the temperature effect in MFP frame- 
work and calculate the Gibbs free energies G of the Hugo- 
niot states for the WZ and RS structures, respectively. 
Fig. 4 shows the Gibbs free energies of these two struc- 
tures as a function of the Hugoniot pressure. It can 
be seen that the overall slope and dependence of Gibbs 
free energy upon pressure for both two structures are 
essentially the same as zero-temperature case [Fig. 3]. 
The Hugoniot WZ^RS transition pressure is 103.8 GPa. 
Compared to zero-temperature result (105 GPa), one can 
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FIG. 4: (color online). Calculated Gibbs free energies of the 
WZ and RS type BeO as a function of the Hugoniot pressure 

see that the temperature effect (typically of 1000 K, see 
Fig. 7 below) is small for BeO WZ— >RS phase transi- 
tion. Thus we derive a conclusion that the enthalpy can 
be reasonably used to study the phase transition of BeO 
in experimentally relevant temperature region. This is 
not surprising considering the fact of high melting point 
of BeO. 

Now we turn to study the Hugoniot for BeO. Hugo- 
niot states, which are derived by the conventional shock- 
wave technique^, are characterized by using measure- 
ments of shock velocity (U s ) and particle velocity (U p ) 
with Vh/Vq = {U s - U p )/U s and P H = PoU s lI p , where 
Ph is the Hugoniot pressure, and po is the initial den- 
sity. By applying the Rankinc- Hugoniot jump conditions, 
these data define a compression curve [volume (Vh) ver- 
sus pressure (Ph)] as a function of known Hugoniot en- 
ergy (E H ): 

Ph(V -V) = 2(E h -E ) (8) 

where Vq and Eq refer to the atomic volume and energy 
under ambient condition, respectively. 

Unlike the static EOS, the temperature along the 
Hugoniot can undergo a range from room temperature 
to several tens of thousands of degrees, thus the cal- 
culations of the Hugoniot state could serve as a good 
check of a theoretical method for the thermodynamic 
calculation. Using the calculated Vq and Eq, we calcu- 
late the Hugoniot EOS for BeO at pressures up to 200 
GPa. Due to the fact revealed in Fig. 4 that there oc- 
curs a WZ— >RS phase transition at the Hugoniot pres- 
sure Pt — 103.8 GPa, thus during calculation of Hugoniot 
EOS, we choose BeO structure to be WZ for the Hugo- 
niot pressure Ph < Pt, while for Ph > Pt the RS phase 
of BeO is chosen. Correspondingly, the initial ambient 
volume (Vq) and energy (Eq) of WZ phase is used when 
Ph < Pt, while for Ph > Pt, the values of Vq and Eq 
refer to the RS phase. Our calculated Hugoniot EOS is 



FIG. 5: (color online). Calculated Hugoniot (solid line) in the 
form of Hugoniot pressure versus Vo/V plot and experimental 
data (squares). 

shown in Fig. 5 (solid curve) in the form of a Hugo- 
niot pressure (Ph) versus Vq/V. Note that to reproduce 
the experimental density^, in accord with the commonly 
used expedient, here the 0-K isotherm of WZ structure 
was uniformly adjusted by reducing 5.85 GPa to each 
pressure and the energies are subsequently modified in 
a consistent fashion. For comparison with our theoreti- 
cal result, the experimental datai^ are also shown in the 
figure (squares). One can see that an overall agreement 
between our calculation and the experiment is fairly good 
up to the experimentally attainable Hugoniot pressure of 
102 GPa, except for four points in the range of 0-30 GPa 
where the strength effects are large. At present, unfortu- 
nately, no experimental shock-wave data at further high 
compression (Ph > 102 GPa) exist. Therefore, to verify 
our theoretical Hugoniot for RS phase at Ph > Pt GPa, 
more experimental data are needed. However, based on 
the good agreement between our calculations and the ex- 
periment for WZ phase BeO, and the coincidence be- 
tween our calculated result of WZ-RS phase transition 
and the previous calculation^, we expect that our calcu- 
lated Hugoniot curve of Ph versus Vq/V can be adopted 
in the realistic application if the higher-pressure data are 
needed. 

For further illustration, we have also transformed the 
calculated Ph versus (Vq/V) relation to that between 
shock velocity U s and particle velocity U p . The result 
is plotted in Fig. 6 (solid curve). The experimental 
data- are also shown in this figure (squares) for com- 
parison. One can see that at mediate particle velocities 
(2 km/s< Up <3 km/s), the calculation agrees well with 
the experiment. Whereas at low particle velocities the 
discrepancy becomes apparent although the calculated 
starting point (U p = 0) coincides exactly with the mea- 
surement (see Fig. 6). Boettger et al. have also noticed 
this discrepancy in their zero-temperature calculation 9 . 
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FIG. 6: (color online). Calculated Hugoniot in the form of U s 
versus U p plot (solid line) and experimental data (squares). 
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FIG. 7: Calculated temperature along the principal Hugoniot 
for BeO. During calculation the crystal phase of BeO has 
been chosen to be wurtzite at pressures lower than WZ-RS 
transition point P t , and rocksalt for P > P t . 



This may be due to the experimental failure to correctly 
measure the particle velocity U p at low shock-pressure 
region by the presence of strength effect as mentioned 
above. 

We also carried out the calculation of the temperature 
along the principal Hugoniot. Note that in the tradi- 
tional reduction of the Hugoniot data, the temperature 
estimate remains less secure since it requires the accurate 
knowledge of specific heat and the Gruneisen parameter 
values that are not well known. Using the present MFP, 
however, all these quantities can be calculated straight- 
forwardly (see below). The calculated Hgoniot empera- 
ture is shown in Fig. 7 as a function of pressure. During 
calculation, again, we choose BeO structure to be WZ at 
Hugoniot pressure Ph < Pt, while for Ph < Pt the RS 
phase of BeO is chosen. The turning segment in Fig. 7 
around Pt is due to this choice of different crystal phases 
in the two pressure regions separated by transition point 
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FIG. 8: Calculated thermodynamic Gruneisen parameter 
along the principal Hugoniot for BeO. During calculation the 
crystal phase of BeO has been chosen to be wurtzite at pres- 
sures lower than WZ-RS transition point Pt, and rocksalt for 
P > Pt- 



Ph = Pt- At present there is no experimental data of 
Hugonoit temperature available. Thus our calculated re- 
sult needs to be verified in the future. 

Since we have explicitly calculated the Hermholtz free 
energy F(V, T) as a function of V and T, all other ther- 
modynamic parameters can be calculated. In particular, 
one can evaluate the thermodynamic Gruneisen param- 
eter. Given the isothermal bulk modulus Bt which is 
given by Bt{V,T) = V ^§y?^ , the volume thermal- 

the constant-volume heat capacity Cy = Cy n {V 1 T) + 
y Ea &r' T ^ ' ^ c thermodynamic Gruneisen parameter 



expansion coefficient (3p(V,T) 



can be obtained as follows: 



7th = 



VB T (V,T)0p{V,T) 
Cv(V,T) 



(9) 



Figure 8 plots our calculated jthVo/V as a function of 
Hugoniot pressure Ph- Again two crystal structures, i.e., 
the wurtzite phase for Ph < Pt and the rocksalt phase 
for Ph < Pt, are used in the calculation. One can see 
that the conventional assumption 7/V=constant for the 
reductions of shock-wave data is not applicable on the 
whole for the present case of BeO. 

Figure 9 plots the constant-pressure heat capacity Cp 
which is calculated by 

C P {V,T) = Cy{V,T) + VTB T (V,T)0 P {V,T). (10) 



For comparison, the constant-volume heat capacity Cy 
and its ionic component Cy n are also shown in Fig. 9. 
One can see that the thermal electronic contribution (the 
difference between the dashed line and dotted line) to the 
heat capacity is negligibly small in a wide range of pres- 
sure. This is what one expects considering the insulator- 
nature of BeO. For the metals, on the other side, the 
thermal electronic contribution to Cy may be compara- 
ble with the lattice ion contribution at high compression 
due to the increasing electronic DOS at fermi level. 
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FIG. 9: (color online). Calculated heat capacities along 
the principal Hugoniot for BeO. During calculation the crys- 
tal phase of BeO has been chosen to be wurtzite at pres- 
sures lower than WZ-RS transition point Pt, and rocksalt 
for P > P t . The solid, dashed, and dotted curves represent 
the calculated constant-pressure specific heat Cp, constant- 
volume specific heat Cv, and the lattice ion only constant- 
volume Cy n , respectively. 



In summary, by using first-principles FPLAPW total- 
energy calculation method supplemented with MFP 
treatment to take into account the vibrational contribu- 
tion of the lattice ions, we have systematically studied 
the Hugoniot EOS of shock-compressed BeO. Our calcu- 
lated Hugoniot shows good agreement with the experi- 
mental data. The other thermodynamic properties, such 
as the Hugoniot temperature, the Griineisen parameter, 
and the heat capacities including the ionic and electronic 
contributions, have also been calculated. In addition, 
the 0-K enthalpy and finite-temperature Gibbs free en- 
ergies of BeO with WZ, ZB, and RS structures have also 
been calculated, which indicates that the WZ— >RS tran- 
sition pressure change by Hugoniot temperature is tiny. 
Therefore, the 0-K calculation is considerably reasonable 
in studying the solid-solid phase transition of BeO. We 
expect our present results shed some light on understand- 
ing the Hugoniot properties of oxides under shock-wave 
compression. 

This work was partially supported by CNSF under 
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1 G. Vidal-Valat, J.P. Vidal, K. Kurki-Suonio, and R. Kurki- 
Suonio, Acta Crystallogr. A 43, 540 (1987). 

2 K.B. Joshi, R. Jain, R.K. Pandya, B.L. Ahuja, and B.K. 
Sharma, J. Chem. Phys. Ill, 163 (1999). 

3 R. M. Hazen and L. W. Finger, J. Appl. Phys. 59, 3728 
(1986); S. P. Marsh, High Pressures 5, 503 (1973); also 
see LASL Shock Hugoniot Data, edited by S. P. March 
(University of California, Berkeley, 1980), p. 241. 

4 A.P. Jephcoat, R.J. Hemley, H.K. Mao, R.E. Cohen, and 
M.J. Mehl, Phys. Rev. B 37, 4727 (1988). 

5 K.J. Chang, S. Froyen and M.L. Cohen, J. Phys. C: Solid 
State Phys. 16, 3475 (1983). 

6 K.J. Chang and M.L. Cohen, Solid State Commun. 50, 487 
(1984). 

7 A. Lichanot, M. Caillet, C. Larrieu, R. Dovesi and C. 
Pisani, Chem. Phys. 164, 383 (1992). 

8 A. Lichanot and M. Rerat, Chem. Phys. Lett. 211, 249 
(1993). 

9 J.C. Boettger and J.M. Wills, Phys. Rev. B 54, 8965 
(1996). 

J. C. Phillips, Bonds and Bands in Semiconductors (Aca- 
demic, New York, 1973); J. C. Phillips, Rev. Mod. Phys. 
42, 317 (1970). 

P. E. Van Camp and V. E. Van Doren, J. Phys. : Condens. 
Matter 8 3385 (1996). 
2 C.-J. Park, S.-G. Lee, Y.-J. Ko, and K. J. Chang, Phys. 
Rev. B 59, 13501 (1999). 



13 Y. Cai, S. Wu, R. Xu, and J. Yu, Phys. Rev. B 73,184104 
(2006). 

14 Y. Mori, T. Ikai, and K. Takarabe, Photon Factory Activ- 
ity Report 20(B), 215 (2003). 

15 Y. Wang, Phys. Rev. B 61 Rll 863 (2000); Y. Wang, D. 
Chen, and X. Zhang, Phys. Rev. Lett. 84, 3220 (2000); 
L. Li and Y. Wang, Phys. Rev. B 63, 245108 (2001); Y. 
Wang, R. Ahuja, and B. Johansson, J. Phys.: Condens. 
Matter 14, 10895 (2002). 

16 E. Wasserman, L. Stixrude, and R.E. Cohen, Phys. Rev. 
B 53, 8296 (1996). 

17 J. S. Dugdale and D. K. C. MacDonald, Phys. Rev. 89, 
832 (1953). 

18 Los Alamos Shock Hugoniot Data, edited by S. P. Marsh 
(University of California, Berkeley, 1980). 

19 T. Jarlborg, E. G. Moroni, and G. Grimvall, Phys. Rev. B 
55, 1288 (1997). 

20 P. Blaha, K. Schwarz, G. Madsen, D. Kvasnicka, and J. 
Luitz, WIEN2K, [http: / / www.wien2k.at71 

21 P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 
77, 3865 (1996). 

22 F. D. Murnaghan, Proc. Natl. Acad. Sci. USA 50, 697 
(1944). 

23 N.N. Sirota, A.M. Kuzmina, and N.S. Orlova, Dokl. Akad. 
Nauk SSSR 314, 856 (1990). 



